Prom Individual to Collective Pinning: Effect of Long-range Elastic Interactions 
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We study the effect of long-range elastic interactions in the dynamical behavior of an elastic 
chain driven quasi-statically in a quenched random pinning potential. This is a generic situation 
occuring in solid friction, crack propagation, wetting front motion, ... In the strong pinning limit, 
the dynamic of the chain is controlled by individual instabilities of each site of the chain. Long- 
range correlations in the displacement field and in the force field develop progressively. The system 
self-organizes to a steady state where the propagation of the instabilities is described by scaling 
laws with characteristic critical exponents. These exponents are numerically estimated through 
the analysis of the spatio-temporal correlation in the activity map. Tuning the exponent a of the 
algebraic decay of the elastic interaction with the distance is shown to give rise to three regimes: 
a Mean-Field (MF) regime valid for a < 1 (very slow decay), a Laplacian (L) regime for a > 3 
(rapid decay of interactions), and an intermediate regime 1 < a < 3 where the critical exponents 
interpolate continuously between the MF and L limit cases. The latter regime is shown to display, 
in the range 1 < a < 2, a mean-field type character only for time-correlations but not for space. 
The effect of the driving mode on the avalanche statistics is also analyzed. 

PACS numbers: 05.40+j, 74.60. Ge, 64.60. Ht, 64.60.LxSelf 



I. INTRODUCTION 

Multistability of elastic media in a pinning potential 
is responsible for the complex dissipative behavior ob- 
served in various physical situations like motion of vor- 
tices in type II superconductors!!], dynamics of a ferro- 
magnetic domain wall driven hy an external magnetic 
field through a random mediumo, Charge Density Waves 
(CDWlu, roughening of crack fronts in fracturea, or solid 
frictiorB. 

In the latter case, the dynamic of the slider can be 
reduced to that of the asperities at the surfaces of the 
solidQ. The role of the elastic body is to mediate the 
interactions between asperities and with the pulling ma- 
chine. It has been shown recently that the competition 
between the elastic restoring force due to the bulk and 
the non linear pinning force, due to the contact between 
asperities of different solids, gives rise to multiple stable 
equilibrium positions when the pinning forces are suf- 
ficiently strong, or when the system is sufficiently large. 
This multistability is responsible for hysteretic behaviour 
of asperities, when they are driven quasi-statically over 
the set of pinning centers belonging to the surface of the 



other solid. Dissipation arises from this hysteretic be- 
haviour and it takes place in the solid, which plays as 
well the role of thermal bath. In the case of spreading 
of a partially wetting liquid on a heterogeneous planeQ, 
the evolution of the contact line depends on the com- 
petition between a "pinning force" due to surface het- 
erogeneities, and an elastic restoring force resulting from 
surface tension. In fact, a distortion of wavelength A of 
the three-phase front modifies the liquid-vapor interface 
over the same distance away from the wall, thus result- 
ing, after integration over the whole liquid- vapor surface, 
in a effective clastic restoring force 

f[x)=-k -\ V^T dx i (!) 

J \x-xi\ a 

with a = 2. The local capillary force is a non local func- 
tion of the entire contact line position. The same effect 
holds for crack propagations, where the stress intensity 
factor exhibits a non local dependence on the front ge- 
ometry, with a = 2 as well. In the case of solid friction, 
the couplings due to the three dimensional elasticity are 
scale invariant and their translationaly invariant part has 
the same form (|l|) , with a = 2D — I for a D-dimensional 
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lattice of asperities coupled via a three-dimensional solid 
(thus a = 1 for D = 1 and a = 3 for D = 2, as can-be 
checked by Fourier transforming the usual relational^) . 
Expression (Q) is also compatible with dipolar interac- 
tions (a = 3) in the case of ferromagnetic domain wall 
(D=2), and we will see that it also takes into account 
the usual Laplacian couplings of vortex lattices, charge 
density waves ... 

The effect of the long range coupling on the dynamics 
of the system has been little studied since each event 
involves a large part of the system and can generally 
not be solved perturbatively. However, it is knownllil to 
have a strong influence on the behavior of the system 
as well as on its stability properties. Here, we focus on 
the behaviour of an elastic line (D=l) with long range 
interactions, driven quasi-statically on a substrate with 
quenched disorder, and in the strong pinning limit. We 
want to show the effect of the long-range elastic couplings 
on the fluctuations accompanying the average behaviour 
of the line. 

By analogy with friction, we will call "asperity" each 
point of discretization of the elastic line. In general, two 
situations may hold: if the couplings are small in com- 
parison to the distortions of the pinning potential, then 
the motion of the chain is controlled by the motion of 
each asperity. This is the strong pinning limit. In this 
case, when the local stiffness is sufficiently weak, the sta- 
ble local equilibrium position of an asperity may disap- 
pear, and the "unstable" asperity advances suddenly to 
the next equilibrium position. In contrast, in the weak 
pinning limit, when the couplings are sufficiently strong, 
the line behaves in a first approximation as a whole. In 
such a case, a separation in domains may occur for a suf- 
ficiently large system sizeu. Thus a coarse-graining of the 
system at large enough scales leads to a simple strong- 
pinning regime. The latter case thus appears to be the 
possible relevant limit in the thermodynamic limit of a 
infinite sized system. 

In this paper, we investigate the effect of the range 
of elastic interactions in the strong pinning case, where 
the fluctuating dynamics of the system results from the 
deterministic propagation of local instabilities. More pre- 
cisely, we focus our attention on the interactions between 
asperities and we consider only the elastic displacement 
field components whose wavelength is greater than the 
distance between "asperities" . This approach is comple- 
mentary to that of considering individual jumps without 
couplings!! We show how the nature of the elastic in- 
teractions allows to interpolate between a "delocalized" 
(or mean field) regime and a "local" regime analogous to 
the Laplacian case. The driving mode is shown to affect 
the velocity of propagation of the instabilities and the 
avalanche statistics. 

In the first part of the paper, we propose a discrete 
model of the dynamics which will be compared to the 
continuous one. Then we study in a systematic way 
the dependence on the decay of elastic interactions, a 
in Eq. (Jl|) , of the propagation of instabilities, of the or- 



ganization of the dynamic and of the kinetic roughening. 
The last section is devoted to a discussion about the uni- 
versality of our model, and on the uniformization due to 
clastic interactions. 

Let us recall briefly that a = 1 for an elastic one- 
dimensional chain with "three-dimensional" couplings, 
but a = — 2 for an elastic string with bendingpelasticityEEl, 
a = 2 for fracture dynamics in infinite sglidsa or wetting 
experiments in a free surface geometryu, and a = in 
propagating capillary fronts due to the fluid flow behind 
the frontEj. 

It is important to note that, for a given physical sys- 
tem, a may be changed experimentally by changing the 
geometry of the apparatus. This is, for example, the case 
of wetting, where a = 2 in a free surface geometry, but 
where the couplings are Laplacian at large scales in a 
Hele-Shaw cell (at vanishing capillary number )t3. 



II. THE MODEL 

Our aim is to single out the effect of interactions in 
the propagation of instabilities in an elastic chain (D=l) 
driven quasi-statically along its transversal direction over 
a rough surface. We consider here a discrete one di- 
mensional elastic chain of size L. The chain of size L 
is discretized into L/d blocks, d represents the distance 
between asperities in the strong pinning limit, and more 
generally the distance between "sites". In the case of 
wetting for example, it will be at least on the order of 
the distance between pinning centers. It is not possible 
to have access to information below d in our model; d is 
the lower cut-off of our model and L is the upper cut-off. 

We call x the current site which is the position of the 
site evaluated in the non deformed geometry. The posi- 
tion of the site in the deformed geometry is x' = x+u(x), 
where u(x) is the displacement field of the chain in the 
laboratory frame. We will refer as well to u(x) as the 
difference in the actual position x' of site x to any homo- 
geneous motion of the undeformed geometry. We con- 
sider only the large wavelength components of the dis- 
placement field (A > d), thus u(x) corresponds to the 
displacement averaged over a cell of size d D . In the fol- 
lowing, we take account only for the displacement in the 
direction of pulling, j, transversal to the chain, and we 
neglect the displacement along the chain. The displace- 
ment field is a scalar and it will be noted u[x), as shown 
in Fig. (I). 

Each site of the chain is subjected to a driving force, 
a random pinning force (from the interface) and a long 
range elastic force describing the coupling to the other 
sites of the chain. The continuous equation of motion 
for this system could thus be written in the overdamped 
limit 
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du(x, t) 
dt 



F ext {t) + r)(x,u(x)) 



(2) 
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J \x-x x \ a 

where F ext stands for the driving force and rj(x, u(x)) for 
the quenched pinning force. The left-hand side of the 
equation contains a phenomenological damping term. 

The long range elastic force has been already intro- 
duced in the previous part of this article. Let us remark 
that it may be described by three fundamental proper- 
ties: (1) it is linear in u{x) (linear elasticity), (2) it is 
scale invariant (despite the lower cut-off d and the up- 
per cut-off L), (3) and it is "translationaly invariant". 
This latter property is due to the fact that we call "cou- 
pling forces" only internal forces. These forces do not 
contribute to the average over all the medium. This is 
compatible with the fact that the elastic energy of the 
chain (in absence of pinning) is translationaly invariant. 
Any uniform displacement of the chain does contribute to 
an average restoring force which has to be compensated 
at equilibrium by a pinning force; but contrary to the 
definition used by C.Caroli et al.U, it does not contribute 
to the elastic couplings. These couplings have an infinite 
range due to their algebraic decay. For the sake of sim- 
plicity, we will refer to a as the range of the interaction. 
As a increases, the interaction tends to be concentrated 
dominantly on the nearest neighbors. On the contrary for 
a small a, the interaction tends to be more evenly dis- 
tributed over the system. We will see later how to adjust 
the external driving force to take account for the average 
dynamic of the system. The quenched pinning force, rj, 
represents the interaction between the heterogeneous sur- 
face and the elastic chain. In case of solid friction, it may 
have various physical origins (Hertzian contact between 
asperities of different surfaces, adhesion ...) but it is al- 
ways non linear. Moreover, a large scale description of 
the contact problem supposes that the basic interaction 
is already the combination of multiple elementary pro- 
cesses, giving rise to discontinuities and multivaluationEl. 
At a sufficiently coarse-grained level, the effective force 
of the interaction will loose its continuity. In order to be 
able to describe the possible fixed point of a renormalized 
description of the surface interactions, we choose to use a 
discontinuous pinning force, described only by its statis- 
tical properties. This force is quenched and random. We 
will restrict ourself in the following to white noise. More- 
over, we anticipate that the fluctuations of the amplitude 
of the pinning force are small when the slider moves for- 
ward, and not sufficiently high to allow backward motion 
of the slider. This is also the case when spatial distortions 
of the pinning forces are strong. In this case, it is still 
possible to use the expression (|J) but with the reserva- 
tion that only the positive part of the r.h.s. is considered. 
Thus, the system is supposed to move only forward and 
the pinning force may be described for each site x only by 
a position coordinate i(x) and a threshold value 77^) (a;) 
as shown in Fig. (l.b). In such a manner, our model 
contains the main features of the dynamics of an elastic 
chain, driven over a pinning surface in the strong pin- 



ning limit. It is also believed to capture the large scale 
behavior of any pinning potential. But, considering de- 
liberately the strong pinning regime, we cannot address 
the interesting question of the cross-over scale from weak 
to strong pinning. 

We are interested in the quasi-static limit of the motion 
of the chain. In this case, the damping is supposed to be 
sufficiently strong to allow each site to reach instanta- 
neously its local equilibrium position. In the quasi-static 
approximation, it is well known that, in the stationary 
regime and in the strong pinning limit, the main con- 
tribution to the global displacement is due to hysteretic 
jumps of each element of the chain, resultingrirom the 
multistability of the site and of the delay ruleO. For a 
given external driving, the initiation of a jump is a local 
instability and thus the dynamics cannot be controlled 
simply by that of the center of mass. To avoid the ex- 
plicit introduction of the dynamic of the system (with 
inertia and viscous damping), and thus the necessity to 
integrate the equation of motion, we assume the follow- 
ing conditions. First, the system evolves from one local 
equilibrium position to the next closest one. Second, we 
monitor the driving force F ext (t) actively so that we al- 
ways maintain the system at the edge of stability. More 
precisely, F ext (t) is adjusted so that only one mode is 
metastable (only one site jumps at each time step). Let 
us emphasize that in contrast to many approaches based 
on continuum equations as Eq. (|^), we do not consider 
a constant force but rather a constant but infinitesimal 
mean velocity. In solid friction, it corresponds to a uni- 
form average imposed displacement v on the top side of 
the solid. In wetting experiments, it will correspond to 
motion at an imposed slow rate of flow. Alternatively, v 
may be deduced from a motion where a force is imposed 
on the top side of the solid (given pressure in wetting 
experiments) . v is related, via shear of the solid, to the 
average position of asperities of the chain. When the dis- 
placement is imposed, the discrete version of our model, 
as studied below, belongs to tha-class of extremal mod- 
els (such as invasion percolatior£3, fractuxe modelsa or 
more recently introduced growth modelsEH). With the 
second assumption, the motion will consist in a series 
of equilibrium points, and the equation of motion does 
not have to be integrated. Under the assumption of the 
over-damped nature of the motion, we will discuss a pos- 
sible extrapolation of the obtained results to different 
driving modes. This discussion will allow us to relate 
the avalanche statistics to the observable intermittency 
of the motion of the chain. 

The above introduced model is discretized to allow 
for numerical simulations. The string is described as a 
one-dimensional, periodic array of sites. The unit length 
along the chain is given by the distance between the sites 
of the chain, and the unit time by the time step St el- 
lapsed between subsequent jumps. Our "time" count in- 
deed the number of jumps (that is the main travelled dis- 
tance). When one site depins from one asperity, it jumps 
by a distance Su, whose maximum value gives the unit 
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length in the direction of propagation. The pinning po- 
tential is described by a set of discrete centers, with a ran- 
dom spatial distribution. We assume that the amplitude 
of each pinning force is identical throughout the medium. 
Alternatively, we could have chosen a uniform distribu- 
tion (on a regular lattice) but with a random magnitude. 
We checked that these two variants (and their combina- 
tion) give the same results, as soon as some randomness 
is introduced. For a periodic array of pinning centers 
with the same amplitude, as could have been guessed, 
the chain motion finally locks on a periodic motion of 
little interest. 

The model runs as follows. We start with a uniform 
flat front, u(x) = 0, for each cell. We suppose that these 
positions correspond to equilibrium positions in absence 
of external driving force (that is for each x, rj(x,0) = 0). 
The external driving force F ext is applied progressively. 
The elastic displacement is negligible, thus only the 
pinning force increases first (r/(x,u(x,t)) = —F ext (t)). 
When the load is sufficiently high, one cell overcomes its 
threshold, and jumps. The maximum load is in this case 

Fext = mm x (-r) (x)). 

The cell x* which overcomes its threshold then jumps to 
the next basin i(x*), thus 

u(x*) — > u(x*) + Su 

During the jump, the external driving force is constant. 
The linearity of the problem allows us to calculate the 
new elastic couplings. The elastic forces along the chain 
are modified to 



G(x*,x*) = 1, thus the stiffness for neighboring sites 
is unity. This sets the scale for forces. The jump size 
Su thus reflects simply the distance between asperities 
on the trackEZl. The dynamic is "extremal" in the sense 
that only the site submitted to the maximum force f{x*) 
advances at each time step. 

Fig.|2| shows examples of the space time distribution 
of "active sites" x* in the medium for various couplings, 
a. Fig.[j] shows the time evolution of the displacement 
front along the chain for the same values of a. It appears 
clearly from these figures that the system organizes in a 
statistically stationary state, and that for large a (local 
couplings), the evolution is spatially inhomogeneous. 

We compared several statistical distributions for Su 
e.g. uniform in the interval [a, 1]. Our results at suf- 
ficiently large times and scales did not display any de- 
pendence on a as long as o ^ 1. As noted above when 
a = 1, the pinning centers are periodically distributed, 
with the same strength, and the chain motion ends up 
being periodic with no interesting features, 

Similar models have already been studiedEl in the case 
a = 2 and in the limit as a oo which is the discrete 
version of the Laplacian oneE3c3 and we will refer to them 
in the following. Similarly, the case a = corresponds 
clearly to a mean-field situation and it admits a simple 
analytical solution whichas quite similar to the situation 
solved by Flyvbjerg et alt 4 ]. 

As far as one is interested only in the effect of the long 
range elastic interactions, it appears convenient to con- 
sider a as a continuous parameter of the model, albeit in 
most physical cases, only integer values can be found. 



f(x)^ f(x)+Su.G(x,x*) 

where G(x, y) results from the interaction kernel dis- 
cussed above. The external driving force is then reduced 
to zero in order to prevent another jump, and the exter- 
nal load may be increased until 

Fext = mm x (-r]^ x )(x) - f(x)) 

in order to allow the next jump. The discrete interac- 
tion kernel is chosen as the simplest expression which 
captures the required periodic boundary condition, and 
follows the power-law decay of Eq. (|l|) : 



G(x,x* 



\x — X* 



G(x*,x*) = ~'E X7 l x *G(x,x*) 



for x ^ x* (3) 
(4) 



In the particular case of a = 2, such an expression (||) 
is the exact result of the summation of interactions for a 
periodic chain. 

This describes one elementary step of the model. We 
start from a uniform flat front and run a long sequence 
of such steps until the system reaches a statistically sta- 
tionary state. Condition Eq. (^) stands for the aver- 
age conservation of forces along the chain. We used 



III. KINETIC ROUGHENING 

Starting with an uncorrelated distribution of forces, 
the system organizes, after a transient, into a highly cor- 
related statistically stationary state. One manifestation 
of these long range correlations is the roughening of the 
displacement front of the chain. The roughness of the 
chain is characterized by the scaling of the correlation 
function 



((u(x + i?) - u($)) 2 )o oc x 2C 
or in Fourier space 

(\u(k)\ 2 ) ocfc- 1 -^ 



(5) 



(6) 



Such a power-law behavior of the power spectrum of the 
front position is shown on Fig. |^ for a = 2. (is called 
the "roughness exponent" of the front. From Fig. (||) we 
determine £ = 0.35 ± 0.02 in good agreement with a pre- 
vious determination of this exponent by Schmittbuhl ct 
al 

As can be seen in Fig. (||), £ depends on a. The larger 
a, the larger the roughness exponent, i.e. the more per- 
sistent the front fluctuation. Note that for large a, the 
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roughness exponent exceeds one, and thus Eq. (||) is in- 
appropriate. Either one should use such a correlation 
function computed on the slope of the front (and mea- 
sures C — 1), or revert to the spectral method for deter- 
mining the roughness exponent. For a = 3, we obtain 
£ w 1.2 again consistent with previous determinationscj. 
The latter case, i.e. for local couplings, is equivalent 
to the Edwards- Wilkinson growth model, however, the 
quenched nature of the noise induces a major change in 
the roughness exponent from £ = 0.5 for an annealed 
noise to 1.2 fop-a quenched noise. 

Fisher et alJiS proposed a renormalization group anal- 
ysis of this problem and obtained 



2a -3D 



(7) 



where D is the space dimension, i.e. in our case D = 1. 
Expression (0) has been obtained by replacing the Lapla- 
cian propagator 1/q 2 by l/q a ~ in the original calcula- 
tion in order to take into account the long-range coupling. 
This theoretical prediction is plotted as a dotted line in 
Fig. [jj. We observe a fairly good agreement for 1 < a < 3. 
However, we observe a marked difference for the Lapla- 
cian case, where the theoretical prediction £p = 1 lies 
outside of the error-bars of our measurement. 

The "roughness exponent" can be measured experi- 
mentally. In the case of wetting on heterogeneous-sur- 
faces in a Hele-Shaw geometry, A. Paterson et al.EJ re- 
ported ( = 0.77. Above a scale equal to the aperture 
of the cell, the problem should be described by the case 
a = 3. However, gravity plays a significant role in this 
problem introducing an additional confining term cutting 
down large scale front fluctuations. It is interesting to 
note that when gravity effects are suppressed (horizontal 
cell) the geometry of the invading fluid is similar to that 
of invasion percolation. The relation between invasion 
percolation and the quenched Edwards- Wilkinson prob- 
lem has been discussed by Roux and HansenEHl, in relation 
with a model of this process proposed by CieplacknJ. It 
may be that the expected scaling is restricted to a rather 
narrow window limited from below by the cell thickness, 
and from above by the capillary length (measuring the 
effect of gravity as compared to the surface tension). In 
case of wetting in a free surface geometry (a = 2), E. 
Rolley et al.o have found £ « 1/3, in very good agree- 
ment with our numerical result (C = O.SS^Contrary to 
the explanation proposed by Robbins et al.ca, using Imry 
and Ma arguments, our explanation takes into account 
the local dynamics of the chain, and may be extended 
by the same way to any range of the couplings. Another 
experimental situation which might be compared to our 
computation is the propagation of an interfacial crack. 
Such a situation between two PMMA plates has been 
studied recently by Schmittbuhl et all They measured 
£ ~ 0.55, a value which is much higher than the above 
mentioned value ( = 0.35 for a = 2. The visco-elastic 
mechanical behavior of PMMA, or the initiation of cracks 
ahead of the front which merge with the front giving rise 



to a tortuous geometry where higher order terms may be 
relevant, jaaay be responsible for this difference. Not also 
that RefB3 is irrelevant in our case because the crack is 
not interfacial there. 



IV. DYNAMICAL ORGANIZATION 

Starting with an uncorrelated distribution of forces, 
the system organizes after a transient into a correlated 
statistically stationary state. In particular, the memory 
of the initial state is lost and the distribution of forces 
f(x) reaches a stationary distribution shown in Fig. (|^). 
The forces have a critical value F max depending on the 
value of a. The distribution is asymetric. The higher the 
value of a, the more peaked is the distribution close to 
F m ax- This distribution depends on the spatial correla- 
tions present along the front and on the distribution of 
force increment following an elementary move. The latter 
aspect can be taken into accounted in a mean-field model, 
expected to hold for a < 1 as shown below. The maxi- 
mum depinning force f{x*) as well (related to the exter- 
nal force by f(x*) = — J^*) (£*) — F ext (t)) has a lower 
critical value corresponding to the maximal external force 
for the whole system, F c . For a uniform constant depin- 
ning threshold for example( r]i( x *)(x*) = rj c = —1), the 
lower critical value of the maximum depinning force is 
(1 — F C /L). This has been drawn in Fig. (^]). F c increases 
with a. 



V. PROPAGATION OF INSTABILITIES 

A. Spatio-temporal map of activity 

As shown in Fig. (H) , the active sites are spatially and 
temporally correlated. These correlations can be trivial 
as when a — 0, or more complex for a > 1. In order to 
analyze these correlations for various values of_a, it is of 
interest to study the probability distributionEj p(r, At) 
of having a distance r between the sites active at time t 
and t + At (i.e. the spatial distribution of the Ai th active 
site). From the numerical simulations, we observe that it 
is possible to describe the entire dependence of p(r, At) 
for a > 1 through a scaling form 



Kr,A i ) = At-V^(_l 7 _) 



(8) 



with a dynamic exponent z which describes the spreading 
of the activity with time over a zone of size £ cx At 1 / 2 . 
For a — > co, z = 2, as for a diffusive system. Fig. (]?]) 
shows a data collapse for At ranging from 1 to 64 and for 
a = 2. The exponent z is determined from this collapse 
and it varies continuously with a. The scaling function 
</> displays the following behavior: 



(x) cx 



for x > 1 
for x < 1 



(9) 
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For distances larger than the active length £, the decay 
of the scaling function is characterized by an exponent 
b which also depends on a. Fig. (^J) summarizes the de- 
pendence of b and zona. 

Let us first propose some arguments that allow us to 
understand the different observed regimes as a function 
of a. 

First, it is possible to link the dynamical roughness ex- 
ponent C to the dynamical exponent z in our model. Let 
us consider a starting point (x D , t a ). After a time At, the 
activity as spread over a distance £(At) around x D . The 
number of moves which have been necessary to cover the 
area between the crack fronts at time i and t D + At scales 
as Ax.Au. The front has a self-affine geometry with a 
roughness exponent £, hence Au oc Ax^ . As a result, 
we get the scaling of the time difference with the extent 
Ax = £ as 



hence 



At = e +c 



z = l + C 



(10) 



(11) 



This scaling relation is accurately obeyed in our numeri- 
cal simulations as long as < £ < 1. For larger values of 
C, the above scaling breaks down, and Au oc Ax, hence, 
the effective exponent z — 2 appears. 

Second, if we make the assumption that in the steady 
state there is only weak spatial distortions in the force 
distribution, then p(r, At = 1) should essentially reflect 
the load sharing rule due to the elastic coupling. This 
argument simply predicts 

b = a (12) 

as found in the numerical simulations for a < 3. An 
analogous relation has been found in the propagation of 
initially localized perturbation in the elastic map, with- 
out pinning^J. 

Third, as a increases, the load redistribution is much 
higher for the nearest neighbors than for the rest of the 
chain. More precisely, when a — 3 the interaction force 



fix) = -k 



(u(x) — u(x')) 



r'|3 



dx' 



contains a singular part proportional to 



fsing(x) OC ( ^ 



d 2 U 



(13) 



plus a regular part. The same holds for any higher value 
of a, up to a tending to infinity where f(x) oc d 2 u/dx 2 . 
Hence, the dynamics of the chain is essentially controlled 
by this singular part, and thus is no longer expected to 
depend on a but rather should be equivalent to the case 
a = 3 or a = oc, i.e. the simple Laplacian (local) ker- 
nel, or the Edwards- Wilkinson equation with quenched 



noise universality class. This regime will be referred to 
as Laplacian, or "L" regime, in the following. 

Let us now return to distribution of distances between 
active sites at a time interval At. We observed the exis- 
tence of a region of extent £ oc At 1 / 2 centered on the ini- 
tiation site where most of the activity is confined, we call 
this region "cluster" — although it is not connected in the 
space-time map. These clusters have a self-affine struc- 
ture. Their scaling is indeed given by Ax oc £ oc At 1 / 2 . 
We note that when r > £, the last occupied site is not 
part of the cluster which originates on the former site. 
The way the activity is distributed in space and time is 
controlled by the statistical distribution of the elemen- 
tary jumps between two consecutive active sites which 
displays a very wide distribution, and the temporal cor- 
relations in those jumps. There are two limit cases which 
have been thoroughly explored, and which may serve as 
guides in the analysis: one case focuses on temporal cor- 
relations with a narrow (say Gaussian) distribution of el- 
ementary jumps, and the absence of a typical time scale 
would induce that the activity can be described as a self- 
affine profile in time, with a roughness exponent 1/z. 
The other limit case corresponds to the absence of tem- 
poral correlations focusing on the power-law distribution 
of elementary jump, and where the space-time map of 
activity can be seen as a Levy walk. We do not know 
of any theoretical attempt to combine these two aspects 
to get a general picture. Due to the nature of the ele- 
mentary jump distribution, it is natural to explore the 
second limit case as a reference, and check its domain 
of validity. Thus we now consider the crude hypothesis 
that temporal correlations can be neglected. The prob- 
ability p(r, t) can thus be written as the convolution of 
p(r, 1) with itself r times. For large r, the distribution of 
p(r, t) will therefore converge to a statistical distribution 
which is stable for addition. If the distribution p(r, 1) has 
a finite second moment, then p(r, t) should converge to 
a Gaussian law as a result of the central limit theorem. 
However, this is never the case, since the above scaling 
implies that this would be obtained for a > 3, and we 
have seen that this case is similar to a = 3. In the in- 
termediate case where 1 < a < 3, we have seen above 



that p(r, t) decays as r 



-b 



~ Q , and hence, the second 



moment of r would diverge in an infinite system size. As 
a result, the central limit theorem does not apply, but 
rather p(r, r) converges to a stable Levy law character- 
ized by a power-law decay as r~ a for all r. Hence, the 
power-law tail will be preserved, for large distances, i.e. 
as long as the correlation length is smaller than the sys- 
tem size, £ < L. This is indeed what is observed for 
a > 1. Still in the case where we neglect time correla- 
tions, we can relate the dynamic exponent z to the large 
distance power-law decay through 



1 = a - 1 



(14) 



We see from Figure (^) that away from the value a = 3, 
where z sw 2, the above relation is poorly satisfied. This 



G 



indicates that temporal correlations becomes more and 
more important as a decreases, and cannot be neglected. 

When a approaches 1, all moments of the distribution 
become controlled by the system size, and hence after a 
few time steps, p(r, At) is smeared out over the entire 
domain size and no more power-law tail survives. The 
scaling given by Eq.(^J) indeed breaks down for a < 1. 
As a natural consequence, time correlations also vanish, 
and we enter a Mean-Field (MF) regime, which is inde- 
pendent of a. 

From these arguments, we arrive at a classification of 
different regimes depending on the value of a: 

• a > 3 the Laplacian regime, or the Edwards- Wilkinson 
regime with a quenched noise, where the kernel is equiv- 
alent to a second derivative. 

• 1 < a < 3 an intermediate regime, where the critical 
exponents which are measured depend continuously on 
a. 

• a < 1 the Mean-Field regime, where spatial correla- 
tions are lost in a few time steps, and where — as in the 
L-regime — the value of a does not influence the critical 
exponents. 

B. Activity recurrence 

In the critical steady state, the activity map is highly 
correlated and exhibits scale invariant features in both 
time and space. We studied previously the spatial dis- 
tribution of activity after a fixed time lag At. We now 
turn to a complementary description, focusing on a sin- 
gle site as a function of time. More precisely, if site i was 
active at time t a , we study the distribution, PFiRSr(t), 
of the time delay t, such that the next move at site i oc- 
curs at time t + 1. Fig. (^|) shows such a distribution for 
a = 1, 1.5,2,3 and oo. The very early time behavior is 
dependent on the distribution of displacement u which is 
chosen. If the displacement is distributed over the inter- 
val [a, b) , for large a, we prevent the recurrence of activity 
immediately after a move. However, this effect lasts only 
for a time much lower than the maximum time T* when 
a > 1. It does not affect the distribution over a large 
interval depending on the system size, and thus can be 
discarded from the analysis. Then for a > 1, we observe 
a power-law decay 

PFIRST(t) OCt~ TFIRST (15) 

where tfirst < 2. This power-law terminates at an up- 
per cut-off which scales as the time necessary for a cluster 
to span the entire system size, T* = L z . Let us now sum- 
marize our observations for various a values. 

• a > 3 Consistently with the previous discussion, the 
value of the exponent tfirst is independent of a and 
amounts to tfirst ~ 1-5 in good agreement with previ- 
ous studies. 

• 2 < a < 3 The distribution drops rapidly (faster than 
any power-law) above the cut-off scale T* . The exponent 



tfirst progressively increases and reaches the value 1 as 
a tends to 2. 

• 1 < a < 2 The exponent tfirst saturates to 1, for 
times up to T*, but in contrast to the previous cases, 
the distribution then reaches a plateau from T* to T**, 
before dropping faster than any power-law. The scaling 
of this second time scale (Fig. |l0|) is identical to that of 
T*, T** oc L z . The relative importance of the 1/t region 
and of the plateau, can be estimated by the ratio T* /T** 
which goes from 1 (no plateau) to (no power-law) as a 
decreases from 2 to 1. 

• a < 1 We enter the Mean-Field regime where pfirst is 
constant for times up to T** which is proportional to the 
system size L. Again this is consistent with the mean- 
field regime previously discussed. 

It is interesting to note that the intermediate regime 
is now split in two cases, with a short and long time be- 
havior differing for 1 < a < 2 rather than a continuous 
evolution with an exponent tfirst going from 1 to 0. 
The fiat plateau regime observed for T* < t < T** and 
1 < a < 2 is similar to the Mean-Field regime for a < 1, 
and thus the transition to the MF regime appears not 
to be as brutal at a = 1 as proposed from the scaling 
of the distance, but rather it turns out to be gradual, 
with a mixed regime displaying non-trivial correlation 
for a "macroscopic" time T*, before reaching a Mean- 
Field behavior at later times. The surprising feature is 
that space and time correlations do not disappear simul- 
taneously. It might be interesting to note that the mean 
distance between active sites separated by a delay At 
also changes for a smaller or larger than 2. For a > 2, it 
follows from the previous spatio-temporal analysis (part 
A), that the mean distance between active sites scales as 
the correlation length £ since b > 2. On the contrary, for 
a < 2, it scales as the system size L. 

Let us consider an active site i at time to. After a time 
t, the number of times, n(t), the site i has been active 
can be estimated as the total number of individual step 
proportional to t divided by the number of sites where 
most of the activity is concentrated, i.e. t 1 / 2 , or 

n^oci 1 - 1 / 2 (16) 

Let us now estimate the same number using the distribu- 
tion of activity recurrence T. From the power-law distri- 
bution of T, we can estimate the maximum among n(t) 
numbers — supposed to be uncorrelated — from 

,oo 1 

/ Pfirst(T) dT oc —— (17) 

JT ma!c (n) n[t) 

or T max (n) oc ^/(wast-I). Tj s hrg the fact that t F irst 
is always smaller than 2, we can compute the mean 
time between activity recurrence (T) oc T^~] C FIRST oc 
n (2-Tj7 /fiSr )/(T F/f?sr -i)^ This same average time is also 

equal to (T) = t/n(t). Equating these two estimates 
gives 

n(t) oc t TF1BST - 1 (18) 
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Comparing Eqs. ([l6|) and ((lj), we arrive at the following 
expression for r first- 

2z-l 

tfirst — (19) 

This estimate is based on the fragile assumption that 
time correlations can be neglected. We have seen above 
that a similar hypothesis lead to z = a — 1, from what 
we deduce 

2a -3 . . 

tfirst — — (20) 

a — 1 

This expression however requires that the distribution 
of T is not governed by the upper cut-off L z , and thus 
tfirst > 1 strictly or a > 2 in the above formula. Fig- 
ure ([n]) shows the measured value of tfirst as a func- 
tion of a, together with the above relation. We observe 
a reasonnable agreement for 2 < a < 3. 

C. Avalanche Dynamic 

Up to now, we have considered an ideal driving where 
the front displacement is controlled by adjusting instan- 
taneously the driving force. Assuming an overdamped 
dynamics, we can reconstruct a different type of driving. 

In order to proceed, it is useful to introduce the notion 
of avalanches or bursts. Let us consider the sequence 
of external force F ext (t). Any force F will split the sig- 
nal into consecutive intervals with alternating F < F ext 
domains (or "obstacles"), and F > F ext intervals (or 
"bursts"). Each burst is thus the front motion which 
would result from a constant force F being imposed to 
the system. ■— ■ 

Avalanches were introduced initially^ in order to un- 
derstand some features of invasion percolationEil. How- 
ever, in RefEJ, some approximations were unjustified or 
unfounded. A complete solution of the avalanche statis- 
tics applied to invasion percolation and other extremal 
models is presented in Paczuski et ao. The key prop- 
erty is that it is possible to relate the avalanches in the 
external forcing to the cluster statistics in some extremal 
models such as invasion percolation. The avalanche size 
distribution at fixed F is a power-law distribution up to a 
maximum avalanche size which diverges as F approaches 
to a critical value F c . This mapping allows to derive ge- 
ometrical information on the activity from a simply ac- 
cessible external signal, the driving force. One amazing 
feature is that the latter is one-dimensional but encode 
a multi-dimensional information. 

In our case, it is unfortunately impossible to establish 
such a direct mapping because of the non-local nature of 
the interaction. Hence, the connectedness of the clusters 
is lost, or equivalently, the clusters to be defined from the 
driving force avalanches do not have a straightforward 
geometrical interpretation. Nevertheless, the maximum 
external force F c encountered in the external driving can 



still be identified with a critical point where the corre- 
lation length diverges. This critical point is a depinning 
transition. For a constant driving F strictly smaller than 
F c , the front will advance only over a finite distance and 
stop when it encounters a pinning force larger than F. 
Similarly, if F > F c , the front will never be stopped. 

It is also possible to use the time evolution of the sys- 
tem to define two distributions of avalanche which do not 
require the precise identification of the critical threshold. 
At each time step, t, one can construct the avalanches 
which correspond precisely to the loading F(t). In fact, 
two such avalanches can be considered, for times ei- 
ther larger or slower than t, the forward and backward 
avalanches respectively as shown in Figure [l2|. Figure |l4| 
shows both cumulative avalanche distributions for our 
model with a > 1. They behave as power laws 

N f , b (t) oc t 1 ^" (21) 

with no upper cut-off apart from the system size. 

MaslovEJ showed that for a variety of models, the ex- 
ponent Tf was super-universal, and equal to 2. Conse- 
quently, forward avalanches do not reveal much infor- 
mation on the model. However, in the extremal models 
studied in Paczuski et alE3, the backward avalanche ex- 
ponent is not super-universal and it does give rise to a 
non trivial critical exponent Tf. 

In our numerical results, the universal result Tf = 2 is 
observed for a > 1 . The backward value is also measured 
to be Tf, = 2 for the same values of a. 

Such power-laws with exponents Tf = T(, = 2 also ap- 
pear in a series of uncorrelated random numbers picked 
from the same distribution, however, we cannot infer 
from this result that the forces are uncorrelated. They 
are indeed strongly correlated as can be seen from the 
power spectrum of the signal F ext (t) (Fig.[l5|) which 
shows a power-law extending over long time intervals (yet 
smaller than the time needed to span the system L z ). If 
we were to interpret this signal as self-affine we would 
estimate from such a graph a self-affine Hurst exponent 
C,p — —0.4, a negative exponent which means that the 
signal is strongly anti-correlated. The same power-law 
appears for all values of a larger than 1. One has to 
be cautious while manipulating Hurst exponent which 
are negative, since most properties relying on the self- 
affine nature of a signal break down in this case. In 
fact for most properties, such functions behave as if the 
Hurst exponent was actually zero (white noise). (Oth- 
erwise, in order to reach reliable conclusions, one may 
consider the time-integrated force signal i.e. the energy 
dissipated during the front motion. More precisely, the 
mean force can be substracted from F before integration 
so that only the fluctuating part of the energy is consid- 
ered. For a > 1, this signal has a self-affine exponent 
£ F + 1 ss 0.6, which can be measured using standard 
tools.) The avalanches are thus expected to behave as if 
the force signal was a white noise, hence r/ = t& = 2. 
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For small exponent a < 1, there is a short time (high fre- 
quency) behavior which seems to follow a similar power- 
law behavior. The plateau appearing for very small times 
is probably due only to numerical accuracy and has not to 
be considered. But at large times, the power spectrum 



of the force crosses over to a 



regime which is the 



signature of a (> = 0.5 self-affine signal, analogous to a 
random walk. The avalanche exponents are consequently 
changed to t/,6 = 1.5 ± 0.05 for a < 1. 

Although the previous analysis of avalanches can in 
principle be accessed experimentally, it requires to be 
able to switch from a displacement-controlled to a force- 
controlled mode, what is in general a rather difficult re- 
quirement. There is another, more natural, way to have 
access to the fluctuating driving force through the in- 
termittent front motion. We have seen that a constant 
force control for finite size system either lead to a stop or 
a constant motion, the transition between the two regime 
corresponding to the maximum F ext . Thus a single con- 
figuration will determine the threshold, and for a finite 
size system, this single point generally is not represen- 
tative of the critical behavior, and it contains a rather 
poor information. A way to circumvent this difficulty is 
to introduce a small but non-zero stiffness e in the control 
(Fig-fill), so that as the front advances, the driving force is 
progressively decreased by a quantity e.(u). This insures 
that the motion will never be unlimited. As soon as the 
front is pinned, then the driving force is slowly increased 
up to the depinning limit. If e is small enough, only the 
forces close to the pinning threshold will be probed. In 
addition, the statistics of the front advance between two 
pinned configuration gives some information on either the 
correlations of F ext , or on its statistical distribution close 
to the critical threshold. 

We performed such an analysis for various values of a. 
We observed that the avalanches have a maximum size 
s* which depends on e as 



with 



k w 0.65 



(22) 



(23) 



for all values of a investigated from 1 to 4. This max- 
imum size can be used to scale the avalanche size s 
(Figjlq) , and thus it allows to account for the distribu- 
tion in an e independent manner 



(24) 



where ip is a power-law distribution for small arguments 
tp(x) oc x~ e for x <C 1, and decays rapidly to zero for 
x 3> 1. And s* r is a normalization factor. 

There are different features which may contribute to 
the statistics of s. We here focus on one aspect of the 
problem namely the statistical distribution of the pinning 
force close to the depinning critical point. Note that the 



results which could be derived from the self-affine struc- 
ture of the signal F ext (t) are not valid due to the asym- 
metry of the distribution of forces, p(F ext ) (Fig.||). 
Let us note 



p(F ext ) = (F c - F ext f = 6 fi 



(25) 



the distribution of external forces near the threshold F c 
After n events, the system has lost an elastic force 6, 

S = ne. 

The maximum force encountered follows 
Balancing the two expressions leads to 



(26) 



hence 

_ 1 + 

K ~2 + f3 

From numerical data, we have ~ 1. It leads to 
k « 0.66 



in good agreement with Eq. (|23j) , although temporal cor- 
relations have not been considered here. 



VI. DISCUSSION 

We have presented in this article a study on the cor- 
relations induced by the dampened motion of a finite 
size elastic chain with long range interactions and driven 
along its transversal direction on a surface with quenched 
disorder. 

We have shown that the system organizes after a tran- 
sient in a stationary state with long-range correlations. 
The memory of the initial state is lost. The final state 
does not depend on the initial state, but its characteris- 
tics depend on a. We have pointed out three regimes: For 
very long range couplings (a < £>), a regime controlled 
by finite size effects and analogous to a Mean-Field one. 
For very short range couplings (a > D + 2), a regime con- 
trolled by small distance singularities of the Green func- 
tion and displaying the same characteristics as Laplacian 
couplings. In between, (D < a < D + 2), an interme- 
diate regime where characteristic exponents b, z and £ 
evolve continuously with a. The transition between the 
three regimes of the activity map is a new one. It is 
summarized in figure [l^. The transition to a Mean Field 
uniform regime does not appear simply for a = 3D/2 as 
been proposed by D.Fisher et al.B by Functional Renor- 
malization Group analysis, but it appears clearly from 
the simulations, that spatial and temporal correlations 
do not disappear simultaneously. The transition to the 







Mean Field regime is characterized by a cross-over time- 
length t c growing with a, above which the uniform distri- 
bution is valid and below which the power law behavior 
with non trivial exponent is observed. Even if the strong 
pinning assumption ensures that the instabilities are lo- 
cal, when t c reaches small times, the behavior of the sys- 
tem becomes "delocalized" . The "delocalized" behavior 
refers to a uniform spreading of activity over the system 
above t c . For a > 2, the system exhibits "collective" 
organization (characterized by long-range correlations) 
but "localized" behavior. For a < 1 the system exhibits 
a uniform "delocalized" behavior analogous to the Mean 
Field one. For 1 < a < 2 the system exhibits collective 
organization with distortions at small scales but "delo- 
calized" behavior at large scales (above t c ). Moreover, 
the time t c goes to infinity in the thermodynamic limit: 
it scales as the activity spreading time L z and does not 
introduce a new scaling. Its origin remains unclear to us. 
An analogous transition appears in the external force for 
a < 1 (see Figurep^j) . 

Modifying the parameters of the system, combination 
of the above behavior can also be found. For a mix of var- 
ious couplings for example, the system exhibits a cross- 
over between a regime controlled by the short range cou- 
pling (at small time and small length) and a regime con- 
trolled by the long range coupling (at large time and large 
length) . The transition between these two regimes is size 
independent and the spatial and the temporal transition 
are related by the dynamical exponent of the shortest 
couplings. 

The results presented above have been obtained with 
a discretized system and periodic boundary conditions. 
We have checked that the same results are obtained, 
far away of the boundaries, without periodic boundary 
conditions. Particularly, the exponents obtained are the 
same in this case. The temporal discretization and the 
breakdown of the translational invariance (because only 
one site is allowed to jump at each time step) are not 
limiting. They are consistent with the "strong pinning 
assumption" . At this stage, the simulations do not allow 
us to justify this assumption more strongly. We know 
howeverEJ that there exists a transition from weak to 
strong pinning which could be used to justify a modeling 
at larger scales. But for very long range (small a) and 
small system size we may encounter situation of weak 
pinning, a case out of reach for our model. We suppose 
that the local curvature of the pinning potential is always 
sufficiently high to maintain the system in the strong pin- 
ning limit. In this limit, the results presented here appear 
to be universal, provided some quenched disorder is in- 
troduced: the choice made for the distribution u and for 
the time step only affects short distances. As has been 
mentioned, they are important in the mean field regime, 
in agreement with ref.EI but they do not affect the position 
of the transition to the Mean Field behavior. When no 
disorder is introduced, the system evolves toward a limit 
cycle of period L. The boundary conditions affect the 
distributions only at the proximity of the upper cut-off 



scale, L or L z , and the displacement fields near the bor- 
ders. When disorder is present, it is possible from our 
numerical results to find intermediate scales where the 
behavior of the system is characterized in space and in 
time by power law distributions. The critical exponents 
of the distributions only depend on a as discussed above. 
More surprisingly, due to the temporal anticorrelation in 
the threshold force, the internal avalanche distribution 
seems to be independent of a and does not reflect the 
correlations of the external force. A better accuracy is 
restored with the help of an elastic external driving. This 
gives rise to external avalanches. Let us mention that ve- 
locity has to be defined car.e£ully in the case of driving 
with constant external forceEj. 

Let us turn to the notion of multistability. For a one 
particle system, the multistability is due to the non linear 
concurrent effect of the elastic force and of the pinning 
force. The condition of multistability can be given by 
a linear stability analysis. The particle is unstable when 
the negative curvature of the pinning potential overcomes 
the elastic stiffness. It leads to hysteretic behaviorQ. The 
energy dissipated at very low velocity is related to the 
area of the hysteretic cycle. For a TV-particle system, 
the instability involves more than one particle, due to 
the coupling between asperities. The size of the domain 
involved in the instability depends on the distortions in 
the curvature of the pinning potential and on the range 
of the elastic interactions. In the "weak pinning limit" , 
the one-particle instability condition is not effective, but 
multistability can occur at larger scale. The minimal size 
giving rise to multistability is usually computed by the 
"Larkin Method"lil. This length is associated to a dy- 
namic by blocks whose existence remains controversial 
in the literature. In the "strong pinning limit", the one 
particle case is effective, but the elastic coupling may 
propagate the instability to an avalanche of larger size. 
We showed here that the system thus organizes into a 
"critical state" . There is no intrinsic length scale in such 
states, but the lowest cut-off and the system size. This 
effect gives rise to "kinetic roughening" in finite systems 
and has been well studied in the case of short range elas- 
tic interactions for various type of pinningjij. It does not 
induce block motion involving more than one asperity. 
An interesting question would be: is the "block dynam- 
ics" of the weak pinning limit analogous to the strong 
pinning dynamics? This question remains open. 



VII. CONCLUSION 

We have shown in this paper that in the strong pin- 
ning limit and for a quasi-static driving, a finite elas- 
tic system driven on a disordered surface exhibits long 
range dynamical correlations. This can be shown by the 
roughening of the surface, in agreement with spreading 
or crack front propagation experiments. It can be stud- 
ied numerically by the statistical analysis of the Activity 
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Map. The range a of the interaction allows to interpo- 
late between a Laplacian dynamical regime (for a > a C 2) 
and a Mean Field dynamical regime (for a < a c i). In our 
one-dimensional simulations, with quenched decorrelated 
disorder, a c i = 1(= D) and a C 2 = 3(= D + 2). The tran- 
sition to the Mean Field regime appears to take place 
progressively from the large times when a < 2. Space 
and time correlations do not disappear simultaneously. 

The internal avalanche size distribution seems to be 
"superuniversal" (independent on a) because it does not 
preserve informations on the temporal correlations in the 
threshold force and exhibits the simple l// 2 behavior 
characteristic of jumps. The dependence on a is recov- 
ered by an elastic external loading. 

This study shows one example of physical systems 
where microscopic disorder leads to large scale structures 
as can be revealed by kinetic roughening (large scale dis- 
placement fluctuations) or avalanche distributions. The 
algebraic decay of the interactions determines the criti- 
cal exponents but other details are unimportant in the 
scaling behavior. The large scale power-law behavior is 
a signature of the self-organization of the system. 

This study was designed to describe a "strong pinning" 
limit. It would be of interest to extend the model to the 
weak pinning case. 
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FIG. 1. (a) Schematic motion of the elastic chain in the force field rj(x,u(x)) between t (black) and t + At (grey). The 
redistribution of the elastic forces via the long-range interactions is not visualized in this picture, (b) Schematic description of 
the pinning force for one site x. The abscisse is now the average position of the line. 

FIG. 2. Activity map showing the location of the active sites as a funtion of time for various couplings a = 1(a), a = 2(b) 
and a = 3(c). The system size is L=1024. 

FIG. 3. Displacement front at different time intervals and for various couplings a = 1(a), a = 2(b) and a = 3(c). The 
system size is L=1024. 

FIG. 4. Average power spectrum of the displacement front position in log-log scale for a a = 2 coupling. The dotted line 
shows a power-law fit which corresponds to the roughness exponent £ = 0.35. The system size is L = 1024. 

FIG. 5. Roughness exponent £ versus coupling range a, measured from the power law fits of the displacement correlation 
function. 

FIG. 6. (a) Statistical distribution n(f) of forces along the chain in the stationary regime, for various couplings a. The scale 
for forces is given by the maximum of forces \G(x* , x*).max(u) |. The pinning threshold is taken arbitrarily as r)i< x *\{x*) = — 1. 
(b) Contribution of the maximum depinning forces f(x*) = 1 — F ext (t) (little o) in the distribution n(f) for a = 1. The dashed 
line stands for a fit of n(f(x*)) near the threshold f(x*) c = 1 — F c /L with /3 = 1. 

FIG. 7. Data collapse of seven different probability distributions p(r, At) for a — 2 and time intervals ranging from 1 to 64. 
The system size is L=512. The redistribution exponent is b = 1.95, and the best data collapse is obtained for z = 1.30. 

FIG. 8. (a) Plot of the b exponent, (b) Plot of the z exponent obtained for different values of a. The dotted and dashed 
lines are the proposed asymptotic behavior discussed in the text. The dotted line (2q/3) on the right take account for the 
relation z = 1 + £ and the expression proposed for £ by Fisher et al. 

FIG. 9. Log-Log plot of the first return probability distribution pFiRsr(t) as a function of the time t for various coupling 
range a = 1, a = 1.5, a = 2, a = 3 and laplacian. The system size is L = 512. 

FIG. 10. Data collapse of the first return probability distribution for various system sizes for a = 1.5. The upper cut-off 
scale as L z as well as the extent of the plateau. 

FIG. 11. Plot of the tfirst exponent obtained for different values of a. Two different exponents appear for 1 < a < 2 
in the Mixed Regime. They correspond to the distribution below and above t c . The dotted lines are the proposed asymptotic 
behavior discussed in the text. 

FIG. 12. External loading F ext (t) as a function of time. From this signal, on can define avalanches either by prescribing a 
fixed value of /, and considering the intervals where / remains below this limit, or by constructing these avalanches starting 
from all forcing values F ext {t) with the statistics issued from the model, and choosing either a forward or backward time 
direction. 



FIG. 13. Instantaneous loading F ext (t) as a function of time and external elastic load of stiffness e. 

FIG. 14. Log- log plot of the distribution of the forward and backward avalanches for a > 1. The straight lines are power-law 
fits of exponents —2. 
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FIG. 15. Spectral amplitude of the time fluctuation of the time integration of the threshold force F ex t{t) for various couplings 
a ranging from to 3. The straight line represents a power-law of exponent —2.2 which implies a roughness exponent of the 
direct time force signal of —0.4 for a > 1. The full line represents a power-law of exponent —4 which implies a roughness 
exponent +0.5. 

FIG. 16. Log-log plot of the distribution of external avalanches for a = 2 and for ten different stiffnesses e ranging from 
10~ 4 to 10~ 3 . The straight line is a power-law fit of exponent —1.09. The upper cut-off scales like s* oc e -0 65 . 

FIG. 17. Schematic transition between the three apparent dynamical regimes. Dotted line separates the classical power 
law regime and a new one regime with a "local" behavior at short time and a "delocalized" behavior at large time. The Mixed 
Regime stands for a regime with Mean-field type character only for time correlations and at large times. 
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